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Abstract 

A transversally driven isotropic ferromagnet being under the influence of a static 
external and an uniaxial internal anisotropy field is studied. We consider the dissipa- 
tive Landau-Lifshitz equation as the fundamental equation of motion and treat it in 
1 + 1 dimensions. The stability of the spatially homogeneous magnetizations against 
inhomogeneous perturbations is analyzed. Subsequently the dynamics above threshold 
is described via amplitude equations and the dependence of their coefficients on the 
physical parameters of the system is determined explicitly. We find soft- and hard- mode 
instabilities, transitions between sub- and supercritical behaviour, various bifurcations 
of higher codimension, and present a series of explicit bifurcation diagrams. The anal- 
ysis of the codimension-2 point where the soft- and hard-mode instabilities coincide 
leads to a system of two coupled Ginzburg-Landau equations. 
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1 Introduction 



By now there is a huge amount of work treating pattern formation in magnetically ordered 
substances under equilibrium conditions, i.e. under the influence of static external as well 
as internal fields. It has led to what seems to be both theoretically and experimentally a 
rather complete picture. The situation is much more precarious if the system is driven far 
from equilibrium by applying strong oscillating magnetic fields. Being hailed as a canonical 
example for pattern formation by some authors (comp. e.g. critical questions have 

been raised by others Q . 

The treatment of ferromagnetic resonance phenomena, which has been initiated by Suhl 
and uses spin waves and their nonlinear interactions, has proved to be an important devel- 
opment in this field. Truncating the number of relevant modes properly many results have 
been obtained which explain most of the observed phenomena pretty successfully. They 
comprise the occurrence of threshold values for parametric instabilities as well as typical 
nonlinear effects including various bifurcation scenarios to low dimensional chaos. 

Even the use of the full (nontruncated) set of spin wave equations is doubtful, however, 
if one discusses pattern formation for large aspect ratios in spatially extended systems, 
because it lacks a systematic derivation. Consequently Elmer Q has proposed another 
approach. Starting from the fundamental Landau-Lifshitz equation of motion describing 
the dynamics of the magnetization field on a mesoscopic level and including various damping 
mechanisms (see below) he derived, for fixed values of the system parameters, amplitude 
equations by using multiple-scale perturbation theory. His main objective was to discuss 
his findings as functions of the strength of the driving field. 

We choose as our starting point the dissipative Landau-Lifshitz equation 

S>,f) = -(Sx% + rSx (Sx%) ). (1) 

containing a Gilbert damping term on the right hand side as well. The latter has been 
justified microscopically in [Q, p|. Disregarding the dipolar interaction as far as it cannot 
be incorporated into a local anisotropy field we study it adopting a rather different point 
of view. Based on a work by Traxler et al. Q, where it has been shown that the homo- 
geneous solutions of equation (|l|) possess a surprisingly rich bifurcation scenario when the 
physical parameters of the system are varied, we extend their work by admitting spatial 
inhomogenities. As in their article we keep the amplitude of the driving field fixed and use 
the static external field together with the detuning from resonance as fundamental param- 
eters. Aiming at obtaining a complete overview in the space of parameters we find sub- 
and supercritical soft-mode and, for the first time for this magnetic system, hard-mode 
instabilities also. The corresponding amplitude equations are derived together with all of 
their coefficients explicitly. 

We then proceed to show that our results imply the existence of some what we deem 
to be highly interesting bifurcations of higher codimension. We just mention here the 
codimension-2 example discussed in section 5.3 which arises when the soft- and hard- 
mode instabilities coincide and which leads to a system of two coupled real and complex 
Ginzburg-Landau equations. 
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However, in our opinion the results given here represent a contribution to earning strongly 
driven ferromagnets a solid place among the classic pattern forming forebears in hydrody- 
namics, nonlinear optics, etc. ||. 



2 Equation of motion 

We start from the dissipative Landau-Lifshitz equation of motion (|l]) and assume that the 
effective field H e tt is composed by external fields, i.e. a static field h z and a circularly 
polarized driving field h± as well as internal fields. The latter are an uniaxial anisotropy 
field aS z and an isotropic exchange field (J > 0). 

( d 2 S \ 

H e ff = (h z + a S z ) e z + J I 1 + h± (coswi e x + sinwt e y ). (2) 

To get rid of the explicit time dependence we first transform to a rotating frame of reference, 
a fact which must not be forgotten when interpreting our results. 

S' x = S x cos(uj t) + S y sin(o; t), S' y = —S x sm(u> t) + S y cos(u t), S' z = S z . (3) 

Observing then that eq. (||) preserves the modulus of the local spin density S(x,t), we 
may put |5[ = 1. Eliminating the superfluous coordinate with the help of a stereographic 
projection 

= S' x + iS' y 

1 + S' Z {> 
we obtain the final form for the equation of motion 



5$ 
~dt 



The parameters 



(i-r) 



(S + i^Z + a® (i_^L) -^(1-^) ( 5 ) 



+ J 



2$ 



1 + |$| 2 




5 = h z — 7T and 7 = — (6) 

1 + r 2 ' 1 + r 2 v ' 

denote the detuning from resonance and the renormalized external frequency, respectively. 

Considering the set of parameters characterizing the effective field it seems reasonable to 
scale them in terms of the anisotropy \a\. Thus a takes on only the values ±1, distinguishing 
between the easy-axis (+1) and easy-plane (—1) case. Furthermore in eq. (||) the exchange 
constant J may be absorbed into the variable x. 

Studying the bifurcation behaviour of this system we may confine ourselves to the case 
h± > and 7 > 0. A negative value of h± merely corresponds to shifting the phase of 
the driving field by ir. Furthermore eq. (||) remains invariant under the transformation 
(<3? — > ^,5 — > —(5,7 — > —7), implying that all bifurcation diagrams are symmetric with 
respect to the origin in the £-7 plane. 
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Our analytical calculations of bifurcation lines, amplitude equations, etc. are performed 
without imposing any restrictions on the parameters T, h±, 5 and 7. For numerical pur- 
poses we will always choose T = 0.1 and keep h± fixed at h± = 0.1. This means that 
we investigate the bifurcation behaviour in the remaining 5-j parameter plane. For the 
codimension-2 bifurcation, to be discussed later on, the variation of the coefficients of the 
amplitude equations will be given as a function of h± also. 



3 Spatially homogeneous fixed point solutions 

We start our systematic investigation of eq. (|5|) by analyzing its spatially homogeneous fixed 
point solutions ( d x &o = dt&o =0). A rather complete bifurcation analysis - including 
codimension-3 bifurcations - of the spatially homogeneous system can be found in Q . Thus 
we confine ourselves to a description of some elementary results which are important in 
the present context. The spatially homogeneous time independent solutions of eq. (||) for 
example are shown to be given by the roots of the algebraic equation: 

( 7 2 + a 2 )S 4 z +2a5 S 3 Z + (S 2 + h\ - a 2 - 7 2 ) S 2 -2a5 S z -5 2 = 0. (7) 

Here S z is the cartesian spin coordinate. From the solution of eq. (Q) S x and S y (and 
therefore $) may be determined explicitly. Depending on the values of the parameters this 
polynomial has either two or four solutions. The corresponding regions in parameter space 
are separated by saddle-node bifurcation lines as is shown in figure 1. As the substitution < Fig.l 
a — > —a and S z — > —S z leaves eq. (0) invariant, this bifurcation does not depend on the 
type of anisotropy chosen. Although the saddle-node bifurcations can be calculated easily 
from eq. (0) we skip this calculation here because an equivalent expression is obtained 
within the framework of our stability analysis in the next section. 



4 Stability analysis 

Let $0 denote a spatially homogeneous fixed point solution of eq. (||) . Proceeding as usual 
in discussing linear stability we add an infinitesimal perturbation 5<&(x,t) 

<f>(x,t) = $ + 6$(x,t). (8) 

Its time evolution is governed in first order by the equation 

<5* = (i - V) [ 9l -d 2 x + ig 2 + (g 3 + ig 4 ) C] 5$, (9) 

where C denotes the operator of complex conjugation, i.e. Cf = f and gi, . . . , 54 are given 

by 

51 = S + h ± Re ($ ) + a 1 ~f ™ ~ 'f o1 (10) 

(1 + \®o\ ) 

92 = j + h ± Im($ ) (11) 

43 = jr^f Re ^ (12) 

54 = (1+iw Im( ^- (13) 
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The plane wave solutions of eq. (^) are written in the form 

5<$> q = cos(qx) (S^ + e xt + <$$_ e**) , (14) 



which accommodates for the presence of the complex conjugation operator. This yields 
after some algebra the following secular equation 

A 2 - A tr(q) + det(q) = 0, (15) 

with 

tr (q) = -2 (g 2 + T(g 1 + q 2 )) (16) 

det(q) = (1+T 2 ) l(g l+q *f + gl-gl-gj\. (17) 

Here tr and det denote the trace and determinant of the two by two ^-dependent coefficient 
matrix which arises when eq. (||) is split into real and imaginary part. As is obvious from 
their definitions {g a }, tr and det depend on the fixed point coordinates 3>o(^> 7> h±, a = ±1) 
as well as on the parameters (5, 7, h±, a = ±1, T) and on the wavenumber q. 

The solution $0 is linearly stable if all possible perturbations 5<& q decay exponentially in 
time, i.e. if ReA±(g) < for all q, where X±(q) denotes the two roots of eq. (|l5|) for some 
fixed wavenumber q. As is described in ||] this yields the conditions 

tr (q) < and det (q) > (18) 

which have to be fulfilled simultaneously for all q. Thus there are two different possibilities 
for these stability conditions to become violated, leading to soft-mode and hard-mode 
instabilities. When both of them are violated simultaneously we have a codimension-2 
bifurcation which requires a separate discussion to be given later on. 



4.1 Soft-mode instability 

A soft-mode instability is characterized by a zero eigenvalue A+(g c ) = which happens if 
det(q c ) vanishes. More precisely it occurs if the minimum value of det(q) becomes zero for 
some finite q = q c , while tr(q) remains negative. Hence we get from eq. (|l~7|) 

d q det\ g=qc = (1 + T 2 ) Aq c ( 91 + q 2 c ) = (19) 

det\ q=qc = (1 + T 2 ) (On + q 2 c f + g\ - g\ - gj) = 0. (20) 
Thus for q c ^ we end up with the bifurcation condition 

52 - al -al = o (21) 

defining a surface in the space of parameters. Its boundary results from —g\ = q 2 = 0, 
which leads to the end points of the bifurcation lines in the 5-7 plane for T and h± fixed. 

The second solution q c = of eq. ([19]) just represents the spatially homogeneous saddle- 
node bifurcation mentioned in section 3. It is described by g\ + g\ — g\ — g\ = 0. The 
corresponding bifurcation lines have been depicted already in figure 1. 



5 



4.2 Hard-mode instability 



In the case of a hard-mode instability the eigenvalues X(q c )± = ±ito(q c ) become purely 
imaginary at the instability point. This occurs if tr(q c ), more precisely its maximum value 
at some q = q c , vanishes, while lo^ = det(q c ) remains positive. From eq. fllBD we get 
immediately 

d q tr\ q=qc = -4 T q c = 0, i.e. q c = (22) 

tr\ Qc=0 = -2 (g 2 + T gi ) = 0. (23) 

Hence the bifurcation manifold is given by eq. ( p3[ ) and its boundary by u\ = det = 0. 

The form of eq. (|l6|) implies that the critical modes are spatially homogeneous. Thus there 
are no primary instabilities creating travelling wave states (uj c ^ 0, q c ^ 0). 

4.3 The bifurcation diagram 

To obtain the soft- and hard-mode as well as saddle-node bifurcation lines in the space 
of the physical parameters we proceed as follows: Keeping T = 0.1, a = ±1, h± = 0.1 
fixed, we solve the two (real) fixed point equations -for convenience we use eq. (|5|), and 
not the polynomial (pi)- together with the proper bifurcation condition numerically. Using 
the PiTCON-package [|io|1 we track the solution curves of this algebraic problem in the four 
dimensional space spanned by Re (<&o); I m (^o)> ^ an d 7- Its projections onto the 5-j plane 
are the bifurcation lines sought for. For the special case of the spatially homogeneous saddle- 
node bifurcation the result has been presented already in figure 1. The soft- and hard-mode 
instabilities which complete this diagram are drawn in figures 2a and 2b for both types < Fig. 2 
of anisotropy. The dotted parts of the bifurcation lines indicate that the corresponding 
instability occurs only after the solution $o being concerned has been destabilized already 
by another type of instability. It is important especially for the hard-mode instability to 
keep this physically "hidden" part of the bifurcation line in mind when discussing some of 
the results which are related to the amplitude equations. 

The gross features^] of the solutions of our system following from figures 2a and 2b may be 
summarized as follows: 

In the case of an easy-axis anisotropy (cf. fig. 2a) we find in region (I) of the parameter 
plane a stable solution near the north pole (3> = 0) and an unstable solution near the 
south pole ($ = oo). Crossing the hard-mode bifurcation line from region (I) to region 
(II) the southern fixed point is stabilized and an unstable spatially homogeneous limit-cycle 
emerges. This limit-cycle in turn destabilizes the northern fixed point solution at the second 
hard-mode bifurcation line. Thus in region (III) we have a stable solution on the southern 
and an unstable solution on the northern hemisphere. Crossing the soft-mode bifurcation 
line separating regions (II) and (IV) the northern (stable) spatially homogeneous fixed 
point solution gets destabilized with regard to spatially periodic states. 

For an easy-plane anisotropy (cf. fig. 2b) the situation in regions (I) and (III) is the same as 
in the easy-axis case. In region (II), however, both homogeneous fixed point solutions are 

1 See reference |7J for more details concerning the spatially homogeneous system. 
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unstable. The northern fixed point of region (I) is now affected by the upper hard-mode 
line but not by the saddle-node bifurcation (76). It looses stability with regard to spatially 
periodic perturbations when passing the soft-mode line separating regions (la) and (VI). 
Note that in accordance with the stability properties of the fixed point involved in the 
pitchfork bifurcation (6) two stable homogeneous states exist in region (la). 

In this survey we have neither mentioned the fact, that on the lower Hopf line a transition 
from super- to subcritical bifurcations occurs implying saddle-node bifurcations of limit- 
cycles nor referred to the global bifurcations of the homogeneous system found in M. 

We stress, however, that our analysis leads via eq. (|i~9| ) to an explicit expression for q c along 
the soft-mode line. Starting from its origin at the homogeneous saddle-node bifurcation 
line q c turns out to increase monotonously from zero to infinity. The corresponding results 
for lu c are contained in j?] already. 

We finally present in figure 3 a number of typical spectra for several fixed points at salient 
positions in the S-j plane. < Fig. 3 



5 Amplitude equations 

After having calculated the possible instabilities and the corresponding bifurcation lines 
in the space of parameters, we go on by analyzing the behaviour of our system in the 
neighbourhood of these marginal stability lines. This will be achieved with the help of 
a perturbation expansion. A natural expansion parameter e measuring the distance from 
threshold results from writing 

5 = 5 c + e 2 5 2 (24) 
7 = 7c + e 2 72- (25) 

(<5 Cl 7c) are the parameters at criticality whereas ($2,72), which may be thought of as the 
components of a unit vector, characterize the direction in which the bifurcation line is 
crossed transversally. 

Let quite generally uo(x,t) denote the critical mode at instability. In order to apply per- 
turbation theory above threshold we put as usual 

$ = <£> + eA(ex,e 2 t)u (x,t) + e 2 r(x,t,e). (26) 

The scaling of the slow space and time coordinates 

X = e x and T = e 2 t (27) 

is motivated by the quadratic dependence of the spectrum on the wavenumber q near 
criticality. The slowly varying amplitude A(X, T) is determined so that the solution ( |26|) 
contains no secular terms. Inserting the ansatz (p6| ) into the equation of motion (||) this 
requirement leads via solvability conditions in lowest non-trivial order of perturbation the- 
ory to an equation governing the time evolution of A. This procedure is well known and 
has been often applied especially to hydrodynamic problems (cf. ref. |ll], |12]| ). For a review 
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see the extensive paper by Cross and Hohenberg 
Kuramoto and Tsuzuki 



as well as the book by Manneville IT 



14 1 used this technique in connection with reaction-diffusion 



equations and derived the coefficients of the amplitude equations for a large class of funda- 
mental equations exhibiting soft- and hard-mode instabilities. Our equation of motion (||) 
does not fit into their scheme because of the quadratic derivative (^§) 2 - It will affect the 
nonlinear coefficient of the soft-mode amplitude equation. 



Concerning magnetic systems this procedure was first applied by Elmer in the context of 
ferromagnetic resonance. His early work is based on a Landau-Lifshitz equation containing 
a dipolar interaction term whereas dissipative effects are taken into account employing the 
Bloch-Bloembergen damping mechanism Jl. More recently he described pattern formation 



in thin magnetic films studying a Landau-Lifshitz equation with Gilbert damping [15]. In 
contrast to our work he is interested mainly in the effects caused by the dipolar interaction. 

The structure of the amplitude equations is determined completely by the nature of the 
instability considered and the symmetries of the underlying system. They are thus univer- 
sal, usually turn out to be real or complex Ginzburg-Landau equations and are found in 
a wide variety of physical systems (see for example the numerous references in Q). The 
individual physics of the system at hand manifests itself in the explicit values of their coef- 
ficients, more precisely in their dependence on the physical parameters. We aim at giving 
a fairly complete survey in this matter for our magnetic problem. 

In the following section we present our results for the amplitude equations for the various 
instabilities and plot the variation of their coefficients along the corresponding bifurcation 
lines. As the necessary calculations in deriving these coefficients explicitly are very lengthy 
they are deferred to appendix A. 

It has been pointed out (see ref. [|[ and the references therein) that the linear coefficients 
of the amplitude equations (cf. eq.(31,32,H)) can be derived alternatively by differenti- 
ating the spectrum of the linearized equation of motion. In our case this amounts to 
-cf. eq. (PJil 



A* = ^~ 



ax 



S-2 



a 



1 #A 

2 dq 2 



(28) 
(29) 



The subscript means that these expressions have to be calculated at the critical values. 
Note that the linear coefficient depends on the direction in which the bifurcation line in 
the (5-7 plane is crossed when entering the unstable region. We will always choose it to be 
normal to that line. The diffusion coefficient a marks the dissipative part of the amplitude 
equation, whereas the nonlinear coefficient r characterizes the type of bifurcation, i.e. for 
example whether it is super- or subcritical, and fixes the saturation value of the amplitude 
in the supercritical case. 
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5.1 Soft-mode instability 



For a soft-mode instability the slowly varying amplitude A(X, T) has been introduced in 
appendix A according to 

$ - $ = e v ( A(e x, e 2 t) e* Qc x + c.c.) . (30) 

Its time evolution is found to be governed by a Ginzburg-Landau equation |8], |ll|, |l2| with 
real coefficients: 

d T A = fi a A + a a d 2 x A-r a \A\ 2 A. (31) 

This equation possesses a Lyapunov potential and is therefore of relaxational type. As the 
modulus of all of its coefficients may be rescaled to unity by appropriate redefinitions of 
time, space and amplitude, the most interesting traits of this equation show up in situations 
where these coefficients change sign and as a consequence the behaviour of the solutions 
changes qualitatively. 

The variation of the coefficients along the bifurcation line is shown in figure 4 for the 
easy axis case. The diffusion coefficient a a , which is proportional to q 2 , is zero at its < Fig. 4 
origin which lies on the saddle-node bifurcation line and increases monotonously. At the 
origin itself the derivation of the amplitude equation looses its validity, an issue which 
will be discussed separately. Notice that the third-order coefficient changes sign along the 
bifurcation line signalling a transition from a subcritical to a supercritical bifurcation. In 
the supercritical case we have a (forward) bifurcation of a stable spatially periodic state, 
whereas for subcritical behaviour the system shows a (backward) bifurcation of an unstable 
spatially periodic state. Whether or not starting from a perturbation theoretical treatment 
in a neighbourhood of this degenerated point higher order terms lead to a stabilization 
remains undecided. In view of appendix A such calculations seem to be impracticable 
however, so that probably one will have to resort the full eq. (||) again to answer such 
questions. Work in this direction will be reported. 

In order to achieve some qualitative understanding of this transition we have performed 
direct numerical simulations on the equation of motion (|l]) using parameter values in the 
neighbourhood of this point. Crossing the soft-mode bifurcation line at point A (cf. fig. 5) < Fig. 5 
in the supercritical region (r > 0) a periodic state emerges continuously. Pursuing that 
solution along the curve C we find that it continues to exist with finite amplitude when 
point B on the subcritical side is reached. So, beyond point B a stable homogeneous and a 
stable periodic solution coexist. We conjecture that -in close analogy to degenerated Hopf 
bifurcations- there is a line of saddle node bifurcations of spatially periodic solutions on 
which the stable periodic solution is destroyed together with the unstable periodic one 
which is created when crossing the soft-mode bifurcation line at point B. 

We have not plotted the coefficients as functions of the parameters for negative anisotropies 
along the segment of the soft-mode bifurcation line connecting the (homogeneous) saddle- 
node and the hard- mode bifurcation lines because this bifurcation is subcritical throughout. 
Recall that all these lines affect one and the same solution. Numerical simulations reveal 
that the solution merely relaxes towards the other stable (spatially homogenous) state after 
a transient period during which periodic structures appear. 
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The obvious divergences of the linear and nonlinear coefficients \x and r are rooted in the fact 
that the critical wavenumber q c vanishes, when the soft-mode bifurcation approaches the 
saddle-node bifurcation. Firstly, the fixed point equations, depending on the parameters S 
and 7, have to be expanded with respect to e 2 in a neighbourhood of the critical parameter 
values. This can be done consistently only if the corresponding Jacobian has non-zero 
eigenvalues. At a saddle-node bifurcation this condition is violated and the expansion 
breaks down leading to a diverging behaviour for the coefficient fi. Secondly, as far as 
the nonlinear coefficient r is concerned, we have to solve an equation second order in e. 
(Refer to appendix A for details). The corresponding solution looses its validity however 
for vanishing critical wavenumber because the linear operator on the left hand side simply 
reduces to the Jacobian so that trivially the solvability condition is no longer satisfied. 
Thus at this point in parameter space the formalism used to derive the amplitude equation 
becomes invalid. A similar situation has been observed by Kuramoto et. al. in ||. 



5.2 Hard-mode instability 

The amplitude equation for this type of instability is given by 

d T B = (fi bR + i m ) B + a b (l + i Cl ) 8 2 x B - r b (l - ic 2 ) \B\ 2 B. (32) 

In contrast to the soft-mode case this equation has complex- valued coefficients. The real 
part a b of the diffusion coefficient is identical to the damping constant T as is proved in 
appendix A, eq. (|98|). The imaginary part \i b j of the linear coefficient can be removed by 
a transformation to a rotating frame [B — > B e lfJ,blt ). 

This complex Ginzburg-Landau equation has been addressed to frequently in the literature 
and is analysed in a wide variety of contexts. The behaviour of its solutions ranges from 
relaxational dynamics for vanishing imaginary parts of the coefficients (c\,C2 = 0) over 
some very complex - even chaotic - types of spatio-temporal patterns in an intermediate 
range of (ci, C2) to the completely integrable limit (c\,C2 — > 00) of the nonlinear Schrodinger 
equation. 

For easy-axis anisotropy the coefficients are plotted as functions of 5 along the whole lower 

hard-mode bifurcation line (cf. figure 2a) in figure 6. The singularities appearing when the < Fig. 6 

bifurcation line approaches the saddle-node bifurcation are caused by the vanishing of the 

critical frequency u> c . Their explanation is completely analogous to that in the soft-mode 

case for vanishing q c in the last section. Regarding the coefficient r b we recognize a change of 

sign signalling a transition from a subcritical to a supercritical bifurcation. Unfortunately 

this degeneracy is found in a region of parameter space, where the corresponding solution 

<3?o is already unstable due to the soft-mode instability separating regions (II) and (IV). 

If the dependence on the driving field h± is taken into account additionally, this situation 

changes as will be expounded in the following section in which the codimension-2 bifurcation 

is treated. 

As has been demonstrated in j7j for the spatially homogeneous system explicitly, and is 



well known under general circumstances [16|, such degeneracy of a hard-mode bifurcation 



entails a saddle-node bifurcation of (spatially homogeneous) limit-cycles. The resulting 
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bifurcation line ends at that point of the hard-mode instability where the coefficient r& 
vanishes. 

Along the upper hard-mode line (cf. fig. 2a) the bifurcation is subcritical everywhere. 
Hence the variation of the coefficients is not shown explicitly. The analysis of the spa- 
tially homogeneous system jij makes evident that this subcritical bifurcation leads to the 
unstable, spatially homogeneous limit-cycle existing in region (II) of parameter space. 

Concerning the easy-plane anisotropy, the variation of the coefficients for both hard-mode 
lines is depicted in figure 7. Along the upper curve (cf. fig. 7a) the bifurcation is supercrit- < Fig. 7 
ical for all values of the parameters whereas along the lower bifurcation line we recognize 
a transition from a sub- to a supercritical bifurcation (cf. fig. 7b). Regarding the super- 
critical portion of that line, the small numerical value of the coefficient r b implies a large 
saturation amplitude of the oscillatory solution. We surmise that this is brought about by 
the stereographic projection technique, which leads to a large modulus of <I>o if the fixed 
point in question lies on the southern hemisphere. It is thus an artificial product. 

5.3 Codimension— 2 bifurcation 

We now focus on the codimension-2 point in parameter space where the hard- and soft-mode 
bifurcations fall together. The spectrum of the linear operator at this point is depicted in 
figure 3 (1) and indicates that the degenerated critical mode is now composed of a spatially 
periodic time-independent and a time-periodic spatially homogeneous part. Supplying both 
of them with slowly varying amplitudes A(X,T) and B(X,T), respectively, we prove in 
appendix A. 3 that their dynamics is governed by two coupled (real and complex) Ginzburg- 
Landau equations 

d T A = v a A + a a d 2 x A-r a \A\ 2 A - s a \B\ 2 A (33) 
d T B = {^ bR + iix hI )B + a b {l + i Cl )d 2 x B -r b {l-ic 2 ) \B\ 2 B 
-{sbR + isbi) \A\ 2 B. 

The linear coefficients as well as the nonlinear coefficients r a , r b and C2 are the same as 

those in the uncoupled equations for the corresponding instabilities. This is obvious from 

continuation arguments. Only the coupling coefficients s a and s b = s b R + i s b i turn out 

to be qualitatively new quantities. To get an impression of their behaviour we study their 

functional dependence of the amplitude h± of the transversal driving field which was fixed 

up to now (cf. fig. 8). We desist from discussing the linear coefficients /i a and [i b because it < Fig. 8 

is inevitable to pay attention to the actual direction in which the unstable region is entered 

in this case (cf. fig. 9). This was unnecessary hitherto, because the direction normal to the 

bifurcation line has been a natural choice. 

For both values ±1 of the anisotropy we recognize that the sign of the coefficient r b changes 
for a certain amplitude h±. Thus for this value, giving rise to a bifurcation of even higher 
codimension, the point of degeneracy of the hard-mode bifurcation, separating its sub- 
and supercritical forms, crosses the codimension-2 point where the hard- and soft-mode 
instabilities coincide. 
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For still higher amplitudes h± - and easy-axis anisotropy (a = +1)- an interesting and 

complex situation arises which is illustrated schematically in figure 9. The degenerate < Fig. 9 

hard-mode bifurcation D enters region (II-III, fig. 2a) so that the stable homogeneous 

(north pole) solution may loose stability not only in a subcritical Hopf bifurcation (left to 

D) but also in a supercritical one (between D and C). 

In a neighbourhood of the codimension-2 point C there are two completely different lim- 
iting cases concerning the behaviour of the system. The first one has to be described (cf. 
region IV) by the real equation for the amplitude of the spatially periodic state, the second 
one (cf. region II) by the complex equation -with all its implications- governing the am- 
plitude of the time-periodic mode. Considering fi a and fib as functions of the angle $ the 
coupled system of equations ( |33l) provides a smooth transition between these dynamically 
radically different situations. There are furthermore saddle-node bifurcations of spatially 
homogeneous limit-cycles on a line emanating from D. We will tackle this problem among 
other topics in a subsequent paper. 

Note that in the case of easy-plane anisotropy the codimension-2 bifurcation always involves 
a subcritical soft-mode bifurcation whereas it is supercritical for easy-axis anisotropy. 

6 Summary 

The first goal of the present paper has been to discuss systematically the instabilities of 
the spatially homogeneous states of a strongly driven ferromagnet against inhomogeneous 
perturbations. Based on the work by Q, in which only homogeneous magnetizations were 
treated, we have determined analytically the stability boundaries of all the fixed point 
solutions in the spatially extended system. Apart from the (spatially homogeneous) saddle- 
node and pitchfork bifurcations we have found hard-mode (u c ^ 0, q c = 0) and soft-mode 
(lo c = 0, q c ^ 0) instabilities as well as codimension-2 situations where both instabilities 
coincide. The corresponding bifurcation lines have been presented in a 5-j parameter plane 
at fixed h± for easy-axis and easy-plane anisotropy. 

We then went on to derive the amplitude equations near those instabilities which ultimately 
lead to pattern formation and have calculated all of their coefficients up to third order 
explicitly. Their functional dependence on the system parameters was shown along each 
bifurcation line which possesses at least one supercritical segment. From this variation the 
locations in parameter space where the real part of the third-order coefficients vanishes have 
been determined for both types of instability. At these points the bifurcation behaviour 
changes from sub- to supercritical. For a Hopf bifurcation such a degenerated situation is 
connected with a saddle-node bifurcation of limit-cycles. The respective bifurcation line in 
parameter space has been calculated already in Q. As explained in section 5.1 we expect 
an analogous behaviour involving saddle-node bifurcations of spatially periodic states in 
the soft-mode case. 

For the interesting and more involved case of the codimension-2 bifurcation we have found 
two coupled Ginzburg-Landau equations, one with real coefficients for the soft-mode am- 
plitude and another one with complex coefficients for the hard-mode amplitude. This set 
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of equations describes a smooth transition between two fundamentally different limiting 
cases. The first one is given by the relaxational dynamics of the real Ginzburg-Landau 
equation for the amplitude of the spatially periodic state. The second one consists of the 
spatio-temporal dynamics described by the complex Ginzburg-Landau equation governing 
the amplitude of the time-periodic mode. In an intermediate range (cf. region (III) of 
figure 9) a competition between both behaviours will occur. A pretty complex situation, 
namely a bifurcation of even higher codimension, arises if the value of h± is chosen in such 
a way that the degenerated Hopf bifurcation falls together with the codimension-2 point. 

We emphasize that it was the explicit though laborious determination of the coefficients 
as functions of the physical parameters which has revealed the existence of these rather 
intriguing situations. As stated in the text there are other ones deserving closer attention. 
We just mention two of them. The first one is the saddle-node bifurcation of spatially 
periodic solutions, which is to be expected in the vicinity of the degenerated soft-mode 
bifurcation. The second one arises at the end points of our bifurcation lines. At these 
points, representing further examples of bifurcations of higher codimension, the perturba- 
tion theoretical ansatz including the scaling of space and time will have to be modified 
substantially. We are working on these problems, hoping that our results will be interest- 
ing not only from a more principal point of view, but -in revealing new aspects of pattern 
formation in magnetically ordered substances- also under realistic experimental conditions. 
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Appendix A: Derivation of the envelope equations 

This appendix contains the calculations of the explicit expressions for the coefficients of 
the amplitude equations starting from the basic equation of motion (|5|). Considering their 
extent they are confined to a reasonable measure. As a general reference for the procedure 
applied see for example ||, || |l3| . For the purpose of this section it is convenient to formally 
rewrite the equation of motion using the substitution <I> = $o + 5& and separate the linear 
and nonlinear contributions as 

f 5® = (t - r) N[6$\. (34) 

N[5&] contains all the nonlinear terms in <5<E> while T which is given by eq. (^) collects the 
linear ones. Both of them depend on the fixed point $o- 

As has been exposed in section || we introduce a small parameter e in order to study the 
system in the neighborhood of the marginal stability line by writing the parameters near 
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their critical values (<5 C) 7e) a s 

5 = 5 c + e 2 5 2 (35) 
7 = 7c + e 2 72- (36) 

The solution of eq. (^) may then be assumed to be of the form 

5$ = e $i + e 2 $ 2 + e 3 * 3 - (37) 

The function <I>i contains the slowly varying amplitude A(X, T) together with the basic 
mode becoming unstable. In the case of a soft-mode instability, for example, it is given by 

$! =vq (A(X,T)e iq °* + cc), (38) 



where vq is a complex number of modulus 1. Considering first the linear part of eq. (34) 
several expansions have to be performed: Because the fixed point coordinates $o depend 
on the parameters 5 and 7 via eq. (0) we have to expand them with respect to e 2 

$0 = $00 + e 2 $02- (39) 
The second order contributions are determined by 



4f_ 
de 2 



^02 + -^r- $02 + -7^0 2 + -7TJ2 = 0, (40) 



e=0 



<9$o 9^0 dtf #7 



where / is an abbreviation for the bracket on the right hand side of eq. (|5|), disregarding 
the spatial derivatives. This leads to 



$00 ($2 + » 72) + $02 (5 C + i 7 C + /i±$oo) + (41) 
$ 02 (l-2|$oo| 2 -|$oo| 4 ) -2*oa*§o 



(1+ l^ool 2 ) 2 



0. 



The coefficients {51, . . . , 34} appearing in the operator f depend on the parameters as well 
as on the fixed point coordinates, which themselves are functions of these parameters. Thus 
the second order quantities in 



9o 



gP+fgg, « = 1,...,4, (42) 



are found to be: 



gf ) = 5 2 + h± Re ($02) - jz ■ f ° , 2V , ($00*02 + $00*02) (43) 

(1 + |$oo| ) 

gf ] = 72 + /i±Im($ 2) (44) 

9 3 ] = h I tT 12^3 Re (" $ oo $02 + $00 ^02) (45) 
(1 + I $00 1 J 

^ = n -u 12^3 Im (-$oo*02 + * 3 ,o*o 2 )- (46) 
(1 + 1*00 1 J 

Finally the introduction of slow space X = ex and time T = e 2 t coordinates leads to the 
substitutions 

d x ^d x + ed x , d t ^d t + e 2 dr. (47) 
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Putting all this together we obtain 



with 



r = r + er 1 + e 2 r 2 , (48) 

f o = d t - (i - T) (g[ 0) -dl + l gf + {gf + z gf ) c) (49) 

fx = 2 (i - r) a, 9 X (50) 

f 2 = 9 T - (i - r) } + z ff f + (<?f + i g® )C). (51) 



Similarly the nonlinear part of eq. (|34|) can be expanded as: 

jV[<5$] = e 2 iV 2 + e 3 iV 3 (52) 

with 

^ = ~ (1 + |^ 00 |2)3 [*00*1*2 + $ o(*l*2 + $1*2) - $00^2] + h|_*l*2 + (54) 

[(2 |$oo| 2 - l)^i^ + ^oo^i + 3 ^00^1^1 -^00^?] + 



(1+ |*oo| 2 ) 4 



4$oo <9*i ^9*i , <9$ 2 \ , n *i-*oo*i /^0*i x2 



+ ^ +2 



1+ [$oo| 2 dx \dX dx J (1+ |$oo| 2 ) 2 V dx 
Inserting the expansions (S^,E3j5^) into eq. (|34|) yields in ascending orders of e: 



r $i = (55) 
f o $ 2 = -f 1 $1 + (t - r) N 2 (56) 
f *3 = -f 2 *i-f 1 $ 2 + (i- T) AT 3 . (57) 

Observe that To, containing a linear and an antilinear part, is a linear operator when 
acting on real- valued functions. In order to solve this hierarchy of equations we notice that 
the solution of the first equation is given by the critical mode itself (neglecting transient 
contributions). This implements a certain definite wavenumber and/or frequency into the 
theory and as a consequence the right hand side of these equations, depending on powers 
of this fundamental solution, can be expanded with respect to harmonics of this mode. In 
this way the solution of the problem may be constructed within a space of strictly periodic 
functions in time and/or space. It is established by applying conventional Fredholm theory 
which necessitates the introduction of an appropriate scalar product. An apt choice proves 

t° De AT 

(/, g):=Re ( ±- [ " dx ^ [ " dt f(x, t)g(x, t)) (58) 



with T c = |jp and A c = ^L. For a soft-mode instability the critical mode is time indepen- 
dent and eq. (58) reduces to an integration over x, whereas for a hard-mode instability we 
have to deal only with a time-integration because there is no (fast) spatial coordinate. 
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In all three cases to be discussed the solvability condition explained below is satisfied 
automatically for eq. ([56"1). The amplitude equations are obtained from the third-order 
eq. (|57]) which is solvable only if 

= -(¥<,, f 2 *i ) - (*o, f i $ 2 ) + (*o, (* - r) N 3 ). (59) 

is the left-null-eigenvector of Tq (which is of course the same as the right-null-eigenvector 
of Ig )• This condition, also known as the Fredholm alternative, assures that the right hand 
side of eq. (|57| ) contains no contribution exciting an eigenvector of Tq with zero eigenvalue. 
The ensuing construction guarantees, that no secular terms occur in the perturbation ex- 



pansion. For further details see M and appendix A. 2 of [13]. 



We now split up our discussion which up to now has been quite general according to the 
different instabilities. 



A.l Soft-mode instability 

From eq. (|55|) the critical mode is given by 

= v (Ae iqcX + c.c), (60) 

(°) +i( a {0) + a {0) ) 
v = 93 [92 + ^^ . (61) 

Vrf ) ) 2 + tf ) +^ ) ) 2 

The left-O-eigenvector of Tq reads 

^0 = w (De iqcX + c.c.) (62) 

(0) (0) _ . (0) 

w = (i-T) 92 +g4 193 (63) 

/ rf ) ) 2 + 0/J ,) +^ ) ) 2 

where the (arbitrarily chosen) normalization 

(w ,v ) = l (64) 

proves to be convenient. The complex parameter D represents an undetermined phase 
factor and is used later on for the extraction of the amplitude equation. 

Inserting 3>i into the second order equation, the latter can be solved because there are no 
resonant terms. We get 

$2 = 2 c 20 |^| 2 + c 22 (A 2 e 2i «° x + c.c.) - v\ — L- d x d x (A e**"* + c.c.) (65) 

92 

with 

- / (my+ ; w+Wy K + " w -'^) <66> 
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C20 
C22 

U2b 



1 

1 

fejL 2 
^0 



{Qc + i gf ) ( n 2a + n 2 b) + (g^ + i gf* ) (n 2a + n 2b ) 



(67) 



2 $oo Vp q 2 
1+ |$oo| 2 ' 



(1+ l^ool 2 ) 3 

4c 



[$ 00 ^ + 2$oo|wo| 2 -$oo^o] 



-3 q 2 + i g 2 ) ( n 2a - n 2h ) + (a^ + i g£> ) (n 2a 
2a 



™2f>J 



(70) 



There is an additional, arbitrary contribution proportional to <I>i which belongs to the 
nullspace of the linear operator Tq. It is not considered any further because it does not 
influence the amplitude equation. 

Putting $i and $2 into the third order eq. (|57|), we finally can exploit the solvability 
condition ([59]). Gathering all terms proportional to D we end up with a Ginzburg-Landau 
equation with real coefficients: 



&r A = n a A + a a d\ A - r a \A\ 2 A. 



(71) 



Analogously the terms multiplying D yield the complex conjugate equation. The coeffi- 
cients are given explicitly by 



2 r 



Ma 



a„ 



i + r 

(0) 
92 

1 + r 2 2 

(0) Qc 
92 



(0) (2) (0) (2) (0) (2) 

92 9 2 -9 3 9 3 ~ 9i g\ 



1 + r 2 



(9 { 3 0) ) 2 + (9 ( 2'+9t") 2 



(0) , (Oh 



Re +of +ig® )(n 3a + n 3b )) 



(72) 
(73) 
(74) 



with 



n3a = h±v (2c 2 o + C22) 
4a 



(75) 



n 3 b 



(1+ |<£oo| 2 ) 3 
6 a 

(1+ l^ool 2 ) 4 
8 $00 «o 9c 



[($00 ^o + $00 ^o) (2 c 20 + c 22 ) + ($00 «o - $00 u o) (2 c 20 + c 22 )] + 



[(2 I $00 1 2 - 1) uo«o + $00 4 + 3 $00 v 2 - $ 4 00 ug] 



1+ |$oo| 2 



C22 + 



2?c «b 



2 „,2 



(1+ |$oo| 2 ) 2 



(w - $ 00 u ) 



(76) 



A. 2 Hard-mode instability 

The derivation of the amplitude equations for the hard-mode instability proceeds along 
the same lines as that for the soft-mode case: the space-periodic critical mode has to be 
replaced merely by a time-periodic one. The main difference between the two instabilities 
is caused by the symmetries of the basic equation of motion. For a soft-mode instability 
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the spatial inversion symmetry leads to a Ginzburg-Landau equation with real coefficients. 
In the hard-mode case complex coefficients are admissible because eq. (||) is not invariant 
under time reversal. 



The linear mode at marginal stability is now given by: 

/ r? 'LLO c t i „ „ \ i n , • ( r> ~i uj c t 



with 





$1 = 




1 


VOa 


Nv 




1 


V0b 


N, 


N v 


= \l 2 



2^ „(o) , , f „(o) 



The left-O-eigenvector of To reads 



(°) , rJ°) „(°) 



c.c. 



where 



w 0a 
N„ 



f - [g? + r 5 f + i (-g^ +r g ™ + (i + 1*) 9 D 



,(°) _i_ r 



^ JO)- 



(0) N * 



2(1+^)^ 
and the following normalization has been chosen: 

{W0a,v 0a ) = (w b,VQb) 
(w a,V 0b ) = (w 0b ,V 0a ) 



(77) 

(78) 
(79) 

(80) 

(81) 

(82) 
(83) 
(84) 



: (85) 
= 1. (86) 

As before the second order eq. (^6|) has no resonant contributions and can be solved: 



$2 = c 2n \B\ 2 + c 2p2 (Be luJct + c.c.)+c 2m 2 i (B e tu)c — c.c 



(87) 



with 



02,, 



C2m2 



C2p2 



i + r 



2 r 



Ujf. 



^(l+iT) n 2n + ( 5 r+^D«2n 



(°) J°h 



2u) c (i- T) n 2p2 + (1 + T 2 ) ( i (i - T) gj 0) re 2m2 + 



r„(°) 4.,- J°)^ 

(.53 + 1 9a ) ™2m2 
-1 



3«g 

2 u; c (z - T) n 2m2 



i + r 2 ) (»(»- r) 5 i 0) n 2p2 + ( 5 r +^D«2p 2 



(°) , „-J )^ 



(89) 



(90) 
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n 2n = h± vl a + h_i_ vl b - j £ ($00 vl a + 2 $ 00 v 0a v 0a - (91) 

(1 + |$oorr 

$00 vl a + $00 4, + 2 $00 V 0b %, " $00 <&) 

4a 

™2m2 = h± v 0a v ob - ($00 v 0a V b + $00 v Qa v ob + (92) 



(1 + l$00 



$00 VQa V0b - $00 V 0a VQb) 

^2 2d — 

n 2p2 = h ± ° a ° b - ($00«0a + 2$ 00 ^0a^0a- (93) 

2 (1 + I $00 1 ) 

$00 via ~ $00 - 2 $00 ^06 ^Ofe + $00 4.) 

In agreement with the soft-mode case we have dropped a term belonging to the kernel of 
To- Using this solution the nonlinear contribution iV 3 of eq. ( |57| ) can be expressed as 

N 3 = \B\ 2 (n 3pl (Be i ^ t + c.c.)+n 3ml i(Be i ^ t -c.c.)) + £ n 3 ,e 4 ^'. (94) 

We refrain from presenting the complete expressions for the coefficients n 3p i and n 3m i 
because they are rather lengthy. They are received by expanding the nonlinear part i¥ 3 in 
terms of the real valued functions (Be tuJct + c.c.) and i(Be luJct — c.c). Using them the 
solvability condition for the third order eq. ( |5^ ) can be evaluated. Again we obtain the 
amplitude equation by gathering all terms in front of E: 

d T B = (fi bR + i m ) B + a b (l + i Cl ) 8 2 x B - r b (l - ic 2 ) \B\ 2 B. (95) 

The coefficients of this complex Ginzburg-Landau equation are: 

(96) 
(97) 
(98) 

(99) 

(100) 

(101) 



/J-bR = 


- 9 ? -r 9 S 2 * 




m = 


1 + r2 ( 9 f»( 9 P»-r 9 f)- 9 f 9 f 


- 9 f 9 f) 


Oib = 


r 




a b ■ c\ = 


(i + r 2 )oi 0) 

~2 [( w 0a, (i - T)n 3m i) + (w ob , (i - 




n = 


r )^3 P i)] 


n-c 2 = 


~2 i( w 0a, (* - r) n 3p i) - (w b, (i - 


r)n 3m i)] . 



A. 3 Codimension-2 point 

In this case the bifurcation conditions ( |2l| ) and (23) for a stationary (u c = 0, q c ^ 0) and 
an oscillatory (u> c ^ 0,q c = 0) instability are satisfied simultaneously. This requirement 
distinguishes a single point (<5 C ,7 C ) in the parameter plane and implies the relations 



lo 2 = (1+T 2 )qi (102) 
9? = -ll (103) 
9 { ? = Tq 2 . (104) 
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The critical mode is now composed of two degenerated modes and consequently one has 
to introduce two amplitudes A(X,T) and B(X,T) to describe their slow modulations, 
respectively: 

$1 = v (A(X,T)e i ^ x + c.c.) + v 0a (B(X,T)e iu ° t +c.c.)+v ab i(B(X,T)e iu ° t -c.c.). (105) 

vq, vo a and vob are defined as previously, compare eq.s (|6li[7^j79|) . The fact that two 
independent degenerated critical modes occur entails that two different left-O-eigenvectors 
enter the formulation of the solvability condition 

* a = w (Ee iq < x + c.c.) (106) 
^Ob = w 0a {De'^ 1 + c.c.)+w Qb i(De i ^ t -c.c.). (107) 

E and D are arbitrary complex parameters. Because most of the following expressions are 
very lengthy, we do not present them in detail but restrict ourselves to just describing how 
to get them. Using a software package for symbolic mathematics it is straightforward to 
obtain them explicitly. 

The notation being necessary to keep book of all the different time and space periodic 
functions and their coefficients is as follows: Introducing a real- valued basic set of functions, 
which consists of fourier modes and their harmonics 

A n+ (x) = A n e inqcX + c.c. (108) 

A n _(x) = i(A n e inqcX -c.c.) (109) 

B n+ {t) = B n e inuJct + c.c. (110) 

B n _(t) = i(B n e inuJct -c.c.), n= 1,2,3,... (Ill) 

all functions may be expanded with respect to them. Of course A = A(X,T) and B = 
B(X,T) still depend on the slow variables. Therefore by inserting in a first step the linear 
mode $i into the second-order eq. (56) its right hand side is built up by the host of functions 

{\A\ 2 , \B\ 2 , A 2+ , B 2+ , B 2 _, A 1+ ■ B 1+ , A 1+ ■ B^}. (112) 

Up to an unimportant term belonging to the nullspace of Tq the solution of this equation 
proves to be 

$2 = c 2A (0, 0) \A\ 2 + c 2B (0, 0) \B\ 2 + c 2A (2, 0)A 2+ + c 2B+ (02) B 2+ + (113) 

P2B-(0,2) 5 2 - + C2AB+0-, 1) Ai + Bi+ + P2AB-(1,1) M+ ^1- + 
P2(l,0) d x d x A 1+ . 

Here the indices (fJ-,v) of the coefficients c 2 ... mark spatial and temporal fourier modes, 
respectively. They are calculated easily by inverting the linear operator on the sub-spaces 
spanned by our basic set of functions. Most of these calculations have been performed 
already in the sections dealing with the soft- and hard-mode instabilities. Only terms made 
up of products of spatial and temporal fourier modes require additional computations. 



Putting $1 and $2 into the remaining third-order eq. ( |57|) N3 becomes 

iV 3 = n 3A+ (l,0)\A\ 2 A 1+ + n 3B+ (l,0)\B\ 2 A 1+ +n 3B+ (0,l)\B\ 2 B 1+ + (114) 
n 3B _(0, 1) \B\ 2 Bx_ + n 3A+ (0, 1) \A\ 2 B 1+ + n 3A _(0, 1) \A\ 2 B x _ + . . . (115) 
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where the dots stand for all those modes which are nonresonant because they are either 
higher order harmonics or constants. Evaluating the solvability condition pertaining to 
eq. (|57]) for each of the two left-O-eigenvectors separately we obtain by collecting all terms 
belonging to E or D, respectively, the following equations governing the time-evolution of 
the amplitudes: 

8 T A = {i a A + a a d 2 x A-r a \A\ 2 A - s a \B\ 2 A (116) 
d T B = (p bR + i m )B + a b (l + ic 1 )d 2 x B-r b (l-ic 2 )\B\ 2 B-(s bR + is bI )\A\ 2 B. 

Hence we have found a Ginzburg-Landau equation with real coefficients for the amplitude 
A being coupled to a complex one for the amplitude B. The coupling coefficients are 

s a = -(w ,(i-T)n 3 B+(l,0)) (117) 
s b R = ~\ [(w 0a ,{i- r> 3 A-(0,l)) + (%,(«- r> 3 A+(0,l))] (118) 

s b i = - [(w 0a ,(i-T)n 3A+ (0,l))-{w 0b ,(i-T)n 3A -(0,l))]- (119) 

All the other first and third-order coefficients are the same as in the soft- and hard-mode 
cases as is clear from continuation arguments. 
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Figure captions 



Fig.l Saddle-node bifurcation line of eq. (||) in the 8-j parameter plane at T = 0.1 and 
= 0.1. The number of spatially homogeneous fixed point solutions is indicated. 

Fig. 2 Bifurcation diagrams including the hard- and soft-mode instabilities of eq. (||) at 
T = 0.1 and h± = 0.1, (a) easy-axis, o = +1, (b) easy-plane, a = —1. The region 
near the homogeneous cusp point is magnified in the insert. 

Fig. 3 Spectra of the linearized operator at various points of the bifurcation diagrams 
(cf.fig. 2). (1) degenerated soft- and hard-mode point, (2) homogeneous cusp point, 
(3) end point of the soft-mode line, (4) homogeneous cusp point, (5) Arnold- Takens- 
Bogdanov point, (6) homogeneous cusp point. 

Fig. 4 Coefficients of the soft-mode amplitude equation for a = +1 on the bifurcation line 13 
of fig. 2a. The coefficient r a changes sign at 5 ~ —0.80, 7 ~ 0.18. 

Fig. 5 Schematic bifurcation diagram for the degenerated soft-mode instability. The labeling 
is explained in the text. 

Fig. 6 Coefficients of the hard- mode amplitude equation for a = +1 on the lower bifurcation 
line of fig. 2a. The coefficients c\ and C2 are shown separately in the insert, rj, changes 
sign at 8 ~ -1.10, 7 ~ 0.036. 

Fig. 7 Coefficients of the hard- mode amplitude equation for a = —1, (a) on the upper hard- 
mode line, (b) on the lower hard-mode line of fig. 2b. The coefficients c\ and C2 are 
shown separately in the insert, changes sign at 5 ~ —1.10, 7 ~ 0.036. The region 
near zero is displayed on a smaller scale (of order 10 -4 ). 

Fig. 8 Coefficients of the codimension-2 amplitude equations (a) for a = +1, (b) for a = — 1. 
The coefficients c\, C2 and a a are shown separately in the insert, changes sign at 
(a) h± ~ 0.90, 5 ~ -2.05, 7 ~ 0.21, (b) h± ~ 1.72, 5 ~ -1.12, 7 ~ 0.31. 

Fig. 9 Schematic bifurcation diagram of the codimension-2 bifurcation and the degenerated 
hard-mode bifurcation at a = +1. The labeling is explained in the text. 
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